A nonperturbative calculation of the electron's magnetic moment 
with truncation extended to two photons 



Sophia S. Chabysheva and John R. Hiller 

Department of Physics 
University of Minnesota- Duluth 
Duluth, Minnesota 55812 
(Dated: October 4, 2010) 

Abstract 

The Pauli-Villars (PV) regularization scheme is applied to a calculation of the dressed-electron 
state and its anomalous magnetic moment in light-front-quantized quantum electrodynamics 
(QED) in Feynman gauge. The regularization is provided by heavy, negative-metric fields added to 
the Lagrangian. The light-front QED Hamiltonian then leads to a well-defined eigenvalue problem 
for the dressed-electron state expressed as a Fock-state expansion. The Fock-state wave functions 
satisfy coupled integral equations that come from this eigenproblem. A finite system of equations 
is obtained by truncation to no more than two photons and no positrons; this extends earlier work 
that was limited to dressing by a single photon. Numerical techniques are applied to solve the 
coupled system and compute the anomalous moment, for which we obtain agreement with exper- 
iment, within numerical errors, but observe a small systematic discrepancy that should be due to 
the absence of electron-positron loops and of three-photon self-energy effects. We also discuss the 
prospects for application of the method to quantum chromodynamics. 

PACS numbers: 12.38.Lg, ll.15.Tk, ll.10.Gh, ll.10.Ef 
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I. INTRODUCTION 



A. Motivation 

High-energy scattering experiments have shown conclusively that the strong nuclear force 
is well described by quantum chromodynamics (QCD), with scattering observables computed 
perturbatively. At longer distance scales, where the properties of hadrons themselves are de- 
termined, the effective couplings are strong and nonlinear, making the derivation of hadronic 
properties from QCD a difficult task. 

Nonperturbative calculations are always difficult, but for a strongly coupled theory such as 
QCD, they are worse. For a weakly coupled theory, one can set aside much of the interaction 
for perturbative treatment and solve only a small core problem nonperturbatively. For QED, 
this core problem is the Coulomb problem; when combined with high-order perturbation 
theory, amazingly accurate results can be obtained for bound states of the theory jlj. In a 
strongly coupled theory one cannot make this separation so easily. 

In the work presented here, the purpose is to explore a nonperturbative method that can 
be used to solve for the bound states of quantum field theories. Although the bound states 
of QCD are of particular interest, the method is not yet mature enough for application to 
QCD. Instead, we will continue with the program developed in the earlier work of Brodsky, 
McCartor, and Hiller and more recently continued by us [9-11], and explore the 



method within QED. This provides an analysis of a gauge theory, which is a critical step 
toward solving a non-Abelian gauge theory, such as QCD. 



B. Fock-state expansions and Pauli Villars regularization 

We will use Fock states as the basis for the expansion of eigenstates. Each bound state 
is an eigenstate of the field-theoretic Hamiltonian, and projections of this eigenproblem 
onto individual Fock states yields coupled equations for the Fock-state wave functions. We 
truncate the expansion to have a calculation of finite size. 

The solution of such equations, in general, requires numerical techniques. The equations 



are converted to a matrix eigenvalue problem by some discretization of the integrals [12 



or by a function expansion for the wave functions [13|. The matrix is usually large and 



not diagonalizable by standard techniques; instead, one or some of the eigenvalues and 



eigenvectors are extracted by the iterative Lanczos process [14], [15j . The eigenvector of the 
matrix yields the wave functions, and from these can be calculated the properties of the 
eigenstate, by considering expectation values of physical observables. 

Although this may seem straightforward, the integrals of the integral equations are not 



finite and must be regulated in some way. We use Pauli-Villars (PV) regularization [16| for 
these ultraviolet divergences. The basic idea is to subtract from each integral a contribution 
of the same form but of a PV particle with a much larger mass. This subtraction will cancel 
the leading large-momentum behavior of the integrand, making the integral less singular. 
To make an integral finite, more than one subtraction may be necessary, due to subleading 
divergences. To arrange the subtractions, we assign the PV particles a negative metric. The 
masses of these PV particles are then the regulators of the redefined theory, and ideally one 
would take the limit of infinite PV masses at the end of the calculation. 

Ordinarily, this method of regularization, being automatically relativistically covariant, 
preserves the original symmetries of the theory. However, it may happen that the negative- 



2 



metric PV particles over-subtract, in the sense that some symmetry is broken by a finite 
amount. In such counterterm is needed, or a positive-metric PV particle can be 

added to restore the symmetry. We have shown in |9j] that a positive-metric PV photon does 
restore the correct chiral limit. 

It is interesting to note that the introduction of negative-metric partners has rece ntly 
been used to define extensions of the Standard Model that solve the hierarchy problem 17 . 
The additional fields provide cancellations that reduce the ultraviolet divergence of the 
bare Higgs mass to only logarithmic. This slowly varying dependence allows the remaining 
cancellations to occur without excessive fine tuning. 

A serious complication in the use of Fock-state expansions and coupled equations is the 
presence of vacuum contributions to the eigenstate. The lack of particle-number conservation 
in quantum field theory means that, in general, even the vacuum can have contributions 
from nonempty Fock states with zero momentum and zero charge. The basis for a massive 
eigenstate will include such vacuum Fock states in products with nonvacuum Fock states, 
since the vacuum contributions do not change the momentum or charge. These vacuum 
contributions destroy the interpretation of the wave functions. In order to have well-defined 



Fock-state expansions and a simple vacuum, we use the light-cone coordinates of Dirac [18 



19|. Light-cone coordinates also have the advantage of separating the internal and external 
momenta of a system. The Fock-state wave functions depend only on the internal momenta. 
The state can then be boosted to any frame without necessitating the recalculation of the 
wave functions. 



C. Numerical methods and Fock-space truncation 

The standard approach to numerical solution of the eigenvalue problem is the method 
originally suggested by Pauli and Brodsky [l2[, discrete light-cone quantization (DLCQ). 
Periodic boundary conditions are applied in a light-cone box of finite size, and the light- 
cone momenta are resolved to a discrete grid. Because this method can be formulated at 
the second-quantized level, it provides for the systematic inclusion of higher Fock sectors. 
DLCQ has been particularly successful for two-dimensional theories, including QCD [20] and 
supersvmmetric Yang-Mills theory [2l|. There was also a very early attempt by Hollenberg 
et al. [22[ to solve four- dimensional QCD. 

Unfortunately, the kernels of the QED integral operators require a very fine DLCQ grid 
if the contributions from heavy PV particles are to be accurately represented. To keep the 
discrete matrix eigenvalue problem small enough, we use instead the discretization devel- 
oped for the analogous problem in Yukawa theory 0, suitably adjusted for the singularities 
encountered in QED. 

An explicit truncation in particle number, the light-cone equivalent of the Tamm-Dancoff 
approximation j23j, can be made. This truncation has significant consequences for the renor- 
malization of the theory in particular the uncancelled divergences discussed below. 



It also impacts comparisons to Feynman perturbation theory [26] , where the truncation elim- 
inates some of the time-ordered graphs that are required to construct a complete Feynman 
graph. Fortunately, numerical tests in Yukawa theory 0, H[ indicate that these difficulties 
can be overcome. The tests show a rapid convergence with respect to particle number. 

To carry out our calculation in QED, three problems must be solved, as discussed in Q. 
We need to respect gauge invariance, interpret new singularities from energy denominators, 
and handle uncancelled divergences. Although PV regularization normally preserves gauge 
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invariance, the flavor-changing interactions chosen for the PV couplings, where emission or 
absorption of a photon can change the flavor of the fermion, do break the invariance at finite 
mass values for the PV fields; we assume that an exact solution exists and has all symmetries 
and that a close approximation can safely break symmetries. The new singularities occur 
because the bare mass of the electron is less than the physical mass and energy denominators 
can be zero; a principal-value prescription is used. These zeros have the appearance of a 
threshold but do not correspond to any available decay. The uncancelled divergences are 
handled (as in the case of Yukawa theory [§]), with the PV masses kept finite and the 
finite-PV-mass error balanced against the truncation error. 

For small PV masses, too much of the negatively normed states are included in the 
eigenstate. For large PV masses, there are truncation errors; the exact eigenstate has large 
projections onto excluded Fock sectors. To see the effect of truncation, consider the form of 
the coupling dependence for the anomalous magnetic moment, which is 



where /xpy is a PV mass scale. The contents of the square brackets are absent in the case 
of truncation. When the large-/xpy limit is taken, this expression becomes zero when there 
is truncation and a nonzero, finite number when there is not. In perturbation theory, the 
order-g 2 terms in the denominator are kept only if the order-g 4 terms are kept in the nu- 
merator, and a finite result is again obtained in the large-/zpv limit. However, the truncated 
nonperturbative calculation includes the order-g 2 terms in the denominator but not the com- 
pensating order-g 4 terms in the numerator. The associated error is minimized by keeping 
/ipv as small as possible, but if too small, the errors associated with keeping the unphysical 
PV Fock states in the basis will be too large. A compromise is to be found for a range of 
intermediate values of /xpy for which physical quantities are independent of /xpy. For QED, 
we see this in the behavior of the anomalous moment of the electron, as a function of the 
PV masses. 

For Yukawa theory, the usefulness of truncating Fock space was checked in a DLCQ 
calculation that included many Fock sectors. The full DLCQ result was compared with 
results for truncations to a few Fock sectors for weak to moderate coupling strengths and 
found to agree quite well We can see in Table 1 of Q that probabilities for higher Fock 
states decrease rapidly. This was also checked at stronger coupling by comparing the two- 
boson and one-boson truncations 0. Figure 14 of [s| shows that contributions to structure 
functions from the three-particle sector are much smaller than those from the two-particle 
sector. 

For QED, there has been no explicit demonstration that truncation in Fock space is a 
good approximation; the two-photon truncation considered here gives the first evidence. 
The usefulness of truncation is expected for general reasons, but a physical argument comes 
from comparing perturbation theory with the Fock-space expansion. Low-order truncations 
in particle number correspond to doing perturbation theory in a to low order, plus keeping 
partial contributions for all orders in a. As long as the theory is regulated so that the 
contributions are finite, the contributions of higher Fock states are expected to be small 
because they are higher order in a. Of course, due to limitations on numerical accuracy, 
we do not expect to be able to compute the anomalous moment as accurately as high-order 



aig 2 [+a 2 c? 4 In /x PV H ] 




1 + b x g 2 + 6 2 fl' 2 ln/ipv H 



perturbative calculations [3, H3] . 
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D. Applications of the method 



The PV regular izat ion method has been considered for QED and applied to a one-photon 
truncation of the dressed electron state 

a sea 

. In Feynman gauge, one PV electron and 



one PV photon were sufficient if the PV electron mass is taken to infinity; otherwise, a 
second PV photon is needed to restore the chiral limit [9[. This choice of regularization has 
the convenient feature of not only cancelling the instantaneous fermion interactions but also 
making the fermion constraint equation explicitly solvable. 

The cancellation of the instantaneous fermion interactions occurs because the individual 
contributions are independent of the fermion mass and have opposite signs. The instanta- 
neous contributions usually provide important infrared cancellations, making their absence 
a possible cause for concern, but these infrared cancellations are instead provided by PV 
contributions. The absence of the instantaneous terms is important for the numerical cal- 
culation, because these terms can greatly increase the computational load, and is significant 
compensation for the addition of the PV fields to the basis. 

Ordinarily, in the light-cone quantization of QED j28|, light-cone gauge (A + = 0) must 
be chosen to make the fermion constraint equation solvable; in Feynman gauge, with one 
PV electron and one or two PV photons, the A + terms cancel from the constraint equa- 
tion. Light-cone gauge was considered in [7J, but the naive choice of three PV electrons for 
regularization was found insufficient; an additional photon and higher derivative countert- 
erms were also needed. The one-photon truncation yielded an anomalous moment within 
14% of the Schwinger term (29[. For the two-photon truncation considered here, the value 
for the anomalous moment should be close to the value obtained perturbatively when the 



Sommerfield-Petermann term |30| is included. However, numerical errors will make this tiny 
correction undetectable, and we will focus on obtaining better agreement with the leading 
Schwinger term of a/27r. 

An extension to a two-boson truncation is interesting as a precursor to work on QCD. 
Unlike the one-boson truncation, where QED and QCD are effectively indistinguishable, 
the two-boson truncation allows three and four-gluon vertices to enter the calculation. A 
nonperturbative calculation, with these nonlinearities included, could capture much of the 
low-energy physics of QCD, perhaps even confinement. 

The approach depends critically on making a Tamm-Dancoff truncation to a finite num- 
ber of constituents. For QCD this is thought to be reasonable because the constituent quark 
model was so successful [III]. Wilson and collaborators [25], 32[ even argued that a light-cone 
Hamiltonian approach can provide an explanation for the quark model's success. The recent 
successes of the AdS/CFT correspondence [33[ in representing the light hadron spectrum of 
QCD also indicates the effectiveness of a truncation; this description of hadrons is equivalent 
to keeping only the lowest valence light-cone Fock state. 

At the very least, the success of the constituent quark model shows that there exists 
an effective description of the bound states of QCD in terms of a few degrees of freedom. 
It is likely that the constituent quarks of the quark model correspond to effective fields, 
the quarks of QCD dressed by gluons and quark-antiquark pairs. From the exact solutions 
obtained using PV regularization in the unphysical equal-mass limit , it is known that 
simple Fock states in light-cone quantization correspond to very complicated states in equal- 
time quantization, and this structure may aid in providing some correspondence to the 
constituent quarks. However, the truncation of the QCD Fock space may need to be large 
enough to include states that provide the dressing of the current quarks, and perhaps a 
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sufficiently relaxed truncation is impractical. As an alternative, the light-front PV method 
could be applied to an effective QCD Lagrangian in terms of the effective fields. Some work 
on developing a description of light-front QCD in terms of effective fields has been done by 
Glazek et al. Q. 



E. Other methods 



A directly related light-front Hamiltonian approach is that of sector-dependent renormal- 
ization [24J , where bare masses and couplings are allowed to depend on the Fock sector. This 
alternative treatment was used by Hiller and Brodsky [35) and more recently by Karmanov, 



Mathiot, and Smirnov [36] . In principle, this approach is roughly equivalent to the approach 
used here; however, the authors of [36| neglect the limitations on the PV masses that come 
from having a finite, real bare coupling, as discussed in [35[ and [llj, and do not make the 
projections necessary to have finite expectation values for particle numbers. 

Our method is complementary to lattice gauge theory 3?}|, which has been studied for 
much longer than nonperturbative light-front methods and has already achieved impressive 
successes in solving QCD. However, it is formulated in Euclidean spacetime and has particu- 
lar difficulty with quantities such as timelike and spacelike form factors, that depend on the 
signature of the Minkowski metric. In contrast, in a Hamiltonian approach with the original 
Minkowski metric, a form factor is readily calculated as a convolution of wave functions. 



A related method is that of the transverse lattice [38|, where light-cone methods are 
used for the longitudinal direction and lattice methods for the transverse. It is, however, a 
Hamiltonian approach which results in wave functions. Applications have been to large- N 
gauge theories and mostly limited to consideration of meson and glueball structure. 

Another approach is that of Dyson-Schwinger equations j39|, which are coupled equa- 
tions for the n-point Euclidean Green's functions of a theory, including the propagators for 
the fundamental fields. Bound states of n constituents appear as poles in the n-particle 
propagator. Solution of the infinite system requires truncation and a model for the highest 
n-point function. Again, as in the lattice approach, there is the limitation to a Euclidean 
formulation. 



F. Outline 



The content of the remainder of the paper is as follows. The structure of the light-front 
Hamiltonian and the Fock-state expansion of the eigenstate are presented in Sec. [TTJ where 
these are used to obtain coupled integral equations for the wave functions. Expressions for 
the normalization and anomalous magnetic moment are also given. Section II I II contains the 
discussion of the solution of the coupled equations, with results for the anomalous moment 
presented in Sec. IIVI A summary of the work is given in Sec. |V] Details of the numerics are 
left to Appendices. 
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II. THE DRESSED-ELECTRON EIGENSTATE 



A. Helicity basis 



For calculations with more than one photon in the Fock space, an helicity basis is con- 
venient. The dependence of the vertex functions on azimuthal angle then becomes simple. 
This will allow us to take advantage of cylindrical symmetry in the integral equations, such 
that the azimuthal angle dependence can be handled analytically. There is then no need to 
discretize the angle in making the numerical approximation. To introduce the helicity basis, 
we define the following annihilation operators for the photon field^ 



1 . . 1 . 

ai± = y|( a ro ± a W > a /(±) = y|( a/1 ^ 10,12 > 



The Hamiltonian can then be rearranged to the form 

i,s P 



(2.1) 



(2.2) 



+J2 d - 



k- 



-a} x {k)ai_ x (k) + aj {x) (k)a H x){k) 



+ J2 I dE d lVli{bi, s (p)hM [ v hs(bg) a l(g-p) 

i,j,l,s,\ 

+ v $s(p,q) a ](\)(Q- p) 

+ &L(p)&i,-»(?) 



U ij,-2s(P, gWixig ~P) + Uij'^aLJq - p) 



+ H.c. 



and the vertex functions become 



V87r 3 (g+-p+)' 
e (p 1 =F ip 2 ){q l ± iq 2 ) + m(m^ 



(2.3) 



p + q + ^fq A 
e q 1 ± ig 2 



p -1 



V$(p,g) 
U+ ]± (p,q) = 0, 



V8T^q + \/Q + ~P A 
e p 1 ± ip 2 



U ii± (p,q) = =F 



rrij{p 1 ± ip 2 ) — rriilq 1 ± ig 2 



p+g+vV" -p 4 
e m^g 4 " — rrijp + 



V8tt^ P + q + Vq + - P + 



1 For details of Fey nman- gauge QED on the light front, particularly the notation, see the discussion in Q 



[»]■ 



and [9(. For a discussion of the residual gauge freedom and the projection onto the physical subspace, see 
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The ai± operators are null, in the sense that [ai±(k), <^ V ±{U)\ — 0; however, we do have 
[a /± (fc),a} /=F (fc')] = -8i V 8{k-U)- 

We will study the state of the electron as an eigenstate of this light-cone Hamiltonian. 
In general, the electron is dressed by photons and electron-positron pairs; however, we 
limit the calculation to photons and truncate the number of photons to two, at most. The 
eigenstate is then expanded in terms of Fock states. In order that the Fock expansion be an 
eigenstate of the light-cone Hamiltonian, the Fock-state wave functions must satisfy coupled 
integral equations. The wave functions are also constrained by normalization of the state. 
The anomalous magnetic moment is then calculated from a spin- flip matrix element. In 
the remainder of this section, we collect the fundamental expressions for the Fock-state 
expansion, the coupled equations for the wave functions, the normalization of the wave 
functions, and the anomalous moment. 



B. Fock-state expansion 

It is convenient to work in a Fock basis where V + and V± are diagonal and the total 
transverse momentum P± is zero. We expand the eigenfunction for the dressed-fermion state 
with total J z = ±| in such a Fock basis as 



1^(3) = E * 6 L=(3|0> + E / d kCg(k)bl s (P - k)al(k)\0) (2.4) 



IJ Sfl 



+ E / dhdhC^ s {hM) 1 b\ s (P- h- h)al(kx)al 



(fe)|o>, 



where we have truncated the expansion to include at most two photons. The Z{ are the 
amplitudes for the bare electron states, with i — for the physical electron and i — 1 for 
the PV electron. The C?£ are the two-body wave functions for Fock states with an electron 
of flavor % and spin component s and a photon of flavor j = 0, 1 or 2 and field component 
fi, expressed as functions of the photon momentum. The upper index of ± refers to the J z 
value of ±| for the eigenstate. Similarly, the C^ kg are the three-body wave functions for 
the states with one electron and two photons, with flavors j and k and field components fi 
and v. 

Careful interpretation of the eigenstate is required to obtain physically meaningful an- 
swers. In particular, there needs to be a physical state with positive norm. We apply the 
same approach as was used in Yukawa theory [6[. A projection onto the physical subspace 
is accomplished by expressing Fock states in terms of positively normed creation operators 
a o^t> a 2At> an d &os an d the null combinations = J2 i \ZCia\ and b\ = + b\ s . The b\ 
particles are annihilated by the generalized electromagnetic current tpj^ip; thus, b\ creates 
unphysical contributions to be dropped, and, by analogy, we also drop contributions created 
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by a^. The projected dressed-fermion state is 

l^(£))ph ys = ^(-l)^&S ± (£)|0) (2.5) 

i 

s/x i=0 j=0,2 fc=j/2 ^ ^ fc 

/• 1 j/2+lfc/2+l ^i+l+m^i/L 



•AW ' i=0 j,fc=0,2 i=j/3 m=fc/2 V 7 ^™ a/1 + 

X ^1-^2)4(^)4,(^)10). 

This projection is to be used to compute the anomalous moment. 

Before using these states, it is important to consider the renormalization of the external 
coupling to the charge 0, 40]. We exclude fermion-antifermion states, and, therefore, there 



is no vacuum polarization. Thus, if the vertex and wave function renormalizations cancel, 
there will be no renormalization of the external coupling. As shown in (9|, this is what 
happens, but only for the plus component of the current. Our calculations of the anomalous 
moment are therefore based on matrix elements of the plus component and do not require 
additional renormalization. 



C. Coupled integral equations 

The bare amplitudes Zi and wave functions and that define the eigenstate must 
satisfy the coupled system of equations that results from the field-theoretic mass-squared 
eigenvalue problem 

M 2 

P-\^(P)) = p+WHP))- (2-6) 

We work in a frame where the total transverse momentum is zero and require that this state 
be an eigenstate of P~ with eigenvalue M 2 /P + . The form of P~ is given in Eq. (12. 2p . The 
wave functions then satisfy the following coupled integral equations: 



[M 2 - m*]* = 




+U£ ± (P-q,P)C£(q)], 



M 



, ml + q 2 ± tf + q\ 



l^-v) y 

[5 s , ±1/2 V£ s (P-q,P) + <WC/4_ s (£-g,P; 

3 

+E(- 1 ) a+v / <ti /r^fW [v£(p-<t-<i>p-<t)<r£tf><i) 

trf J ~ v 1 + hio^ - - - 



(2- 



y/1 + 5 bl 5^ 

+UZ(E -q'-q,P- q')cX_ s ( q ', q 
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(1 - 1/1 - 1/2 

V 7 ! + Sji6^ 



yi V2 

E(- X ) a { v 7 ^ fe s (£ g 2 ) c £(g 2 ) 



(2.9) 



+ ve 



A diagrammatic representation is given in Fig. [TJ The first of these equations couples the 





FIG. 1. Diagrammatic representation of the coupled equations (|2.7p . (j2.8|) . and (|2.9p of the text. 
The filled circles and ovals represent wave functions for Fock states; the solid lines represent 
fermions; and the dashed lines represent photons. The crosses on lines represent the light-cone 
kinetic energy contributions, which are summed over all particles in the Fock state. 

bare amplitudes Zi to the two-body wave functions, C^ s . The second couples the C^ s to the 
Z{ and to the three-body wave functions, C^X. S . The third is truncated, with no four-body 
terms, and simply couples Cj£^ to C^ s algebraically. From the structure of the equations, 
one can show that the two-body wave functions for the J z = — 1/2 eigenstate are related to 
the J z = +1/2 wave functions by 

cg; = -c£*, c^ = c^. (2.io) 

This will be useful in computing the spin-flip matrix element needed for the anomalous 
moment. 
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D. Normalization and anomalous moment 



The projected Fock expansion (I2.5P is normalized according to 
In terms of the wave functions, this becomes 

1 j/2+l 



(2.11) 



iH£(-i)^r + £/^£e. 

i s/t J j=Q,2 



1 \i+k 

E E ^C-£(i) 

i=0 fc=j/2 V?fc 



(2.12) 



sfii/ j,k=0,2 



i=0 l=j/2 m=k/2 



1 + Si m S, 



lm,v fit/ 



Using the coupled equations, we can express all the wave functions and C^ ks and the 
amplitude Z\ through the bare-electron amplitude z$. The normalization condition then 
determines z§. For the two-photon truncation, where the wave functions are computed nu- 
merically, the integrals for the normalization must also be done numerically, using quadrature 
schemes discussed in Appendix [A] 

The anomalous moment a» can be computed from the spin-flip matrix element of the 
electromagnetic current J + [41] 



Qx iQy 

2M 



HQ 2 ) = ±^ ± (P + Q)\^\^(P)) 



phys, 



(2.13) 



where Q is the momentum of the absorbed photon, Fi is the Pauli form factor, and we work 
in a frame where Q + is zero. At zero momentum transfer, we have a e = ^(0) and 



, / 1 i/2+1 

E/^Eo EE 

sfi J j=0,2 \i'=0 k'=j/2 

1 j/2+1 , 



:-i) 



xy 







+ i- 



.7=0,2 

d 



(2.14) 



dk T dk 



EE 

i=0 k=j/2 



l\i+k 



C'iks (^) 



■■■ j,k=0,2 



i>=0 l'=j/2m'=k/2 



X 



E 



Va 



d . d 
+ 1- 



dk n ,x dk 



t= 0l =j/ 2m=k/2 V 1 + ^ 



J flU 



In general, these integrals must also be computed numerically. 

The terms that depend on the three-body wave functions are higher order in a 

than the leading two-body terms. This is because f !2.9j) determines C^ s as being of order 
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yfa or e times the two-body wave functions, the vertex functions being proportional to the 
coupling, e. Given the numerical errors in the leading terms, these three-body contributions 
are not significant and are not evaluated. The important three-body contributions come 
from the couplings of the three-body wave functions that will enter the calculation of the 
two-body wave functions. 



III. SOLUTION OF THE EQUATIONS 



A. Integral equations for two-body wave functions 



The first and third equations of the coupled system, (12. 7\\ and (12.91) . can be solved for 
the bare-electron amplitudes and one-electron/two-photon wave functions, respectively, in 
terms of the one-electron/one-photon wave functions. From (12 .7D . we have 



M 2 



+ P+U£ ± (P-q,P)C£(q)], 



3,1,1* 



(3.1) 



and from (I2.9P we have 



y/1 + 5 3 l^ v 



M 2 - 



X 



(l-l/i -l/a) 



2' 



+ V6 



+K,- s (p-q 1 -g 2 ,p-g 2 )c:t- s {% 

v» ls {P-g 1 -g 2 ,P-g 1 )c»%{g 1 ) 



-g x ~g^P- qJCg^qj] } 



(3.2) 



Substitution of these solutions into the second integral equation (12.81) yields a reduced inte- 
gral eigenvalue problem in the one-electron/one-photon sector. 

To isolate the dependence on the azimuthal angles, we use q\ ± iqf = q i ±e ±l<l ' i and qf = 
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ViP + , and write the vertex functions (12.31) as 

v , p _ _ p _ v __l e (gue^ 1 '^) + g 2 ±)g 2 ± + m 8 m tt 

T/(±)/'D D ^ e±i4>2 e ^2± 

^±(£-5 1 -i 2 .£-« 2 ) = o. 



e m a (gi_Le ±i(<?il ^ + g 2 ±) - rn t Q2± 



U- ± (E gl q_ 2 ,E g 2 )-± (p+)5/2 ^_3_ { i-y 2 ) { i- yi - y2 ) 
^il ) (£-9 1 -£ 2 ,£-9 2 ) = 0, 

T j{±)(p_ _ P _ \ = -,- 1 e ^(1 -y 2 ) -m a (l -i/x -y 2 ) 

mTl - - 1 - 2 '- - sJ V + ) 3/2 Vs^ (i-wKi-yi-w) 



The angular dependence of the wave functions is determined by the sum of J z contribu- 
tions for each Fock state. For example, in the case of C\^_ , the photon is created by 
a ](+) = v^ ? 1 _ ^ a i 2 ^' wn ^ cn contributes J z = — 1 to the state, and the constituent electron 
contributes J z = — |; therefore, to have a total J z of +|, the wave function must contribute 
J z = 2, which corresponds to a factor of e 2l< ^ . For the full set of J z = +| wave functions, we 
find 

C±X(q) = VP+C±+(y, q ± ), Cr+(q) = — ^C$.(y, q±), (3.4) 



C±+_{q) = y/P+e»C%±(y, U ), Cr+(q) = -?L-C£;(y, q ± ), (3.5) 

p 2i<j> i 

C^l\q) = -^C^ + (y,q ± ), C^ + (q) = —tfH+fo, q x ). 



The wave functions have different dependence on longitudinal momenta, resulting in different 
powers of P + , which have been explicitly factored out; they cancel against other P + factors 
in the final integral equations. 

The energy denominator of the three-body wave function can be written as 



M 2 _ + (gi± + gk) 2 _ A + l!± _ + 4± (3 g) 

(1 - Vl ~ V2) V\ V2 

M 2 _ m f + gL + gzi + 2 gi£g2± cos(0i - <j> 2 ) _ $ + g 2 ± _ + ql ± 

(i - yi - 3/2) yi 2/2 
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The light-cone volume element dq' becomes ^P + dy'd(p'dq'l. All the angular dependence can 
then be gathered into integrals of the form 



1„ 



2tt ill 



2tt D ajb (q 1± , q 2 ±) + F(q 1± , q 2 ±) cos(0 - <f)') ' 



with In | = 0, 1, 2, 3 and and F defined as 



F>ajb(q±,q'±) 



m 2 a + ql + q'l , $ + q\ , \i\ + qf 



2q±q'± 



+ 



y 



y' 



i-y-y' 

The integral equations for the two-body wave functions then take the form 



2 m^ + q 2 ± tf + qV 



M z - 



1-2/ 



i'j's'v ^° 



1-2/ 



dy'dq'lAfsZ'Ay, qr, y', 4±)C$}Av'> 



(3.7) 



(3.8) 



There is a total of 48 coupled equations, with i = 0, 1; j = 0, 1,2; s = ±~; and /i = ±, (±). 
A diagrammatic representation is given in Fig. [2j 






FIG. 2. Diagrammatic representation of Eq. (|3.9f) of the text. The conventions for the diagrams 
are the same as in Fig. [TJ 
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The first term on the right-hand side of (13.91) is the self-energy contribution 11 



liu>(y, ?±) = £(-i) </+a+ *6 



i -, r\ ni; -\-Ta-i m.-ri> 

1 dx d 2 k ± mimv - 2^3^m a + 



mi+k 2 



X 71 



(3.10) 



l-x 



with 



+ (1 - y)M 2 - t*±3±. 

y 



(3.11) 



The kernels and in the second and third terms correspond to interactions with zero 
or two photons in intermediate states. The zero-photon kernel factorizes as 



J ljl7j's'(y^ y', q'±) = Yl v i%s(yi q± 



(0)M/ 



(-!)■ 



M 2 



m 



(3.12) 



with 



ija+ 



(0)- 



L ' y/y' Vh 



mmr, 



T/(0)(+) 
) v ija+ 



o v. ( . 0)(_) 



[i-y)Vv 



, (3.13) 



and 



o, v.^- 



m a q ± T/ (o)(+) _ n xr(0)(-) 

^3 H „.\ nr.' 1 ija— u ' i'ja- 



')--V)y/V 



^ m f-^- m \ (3.14) 



[i-y)Vy 



The two-photon kernels are considerably more involved, large in number, and not particu- 
larly illuminating; details will not be given here but can be found in 10|. The associated 
angular integrals X n are given in detail in Appendix [Bj 



B. Fermion flavor mixing 

The presence of the flavor- changing self-energies, the I mi with i ^ i', leads naturally to a 
fermion flavor mixing of the two-body wave functions [HI]. The integral equations for these 
functions have the structure 

AojCojs ~~ BjC\fs = ~ QZ^ojs^ (3.15) 

where and Bj are defined by 

Aij = m? + gi + + ±J&- - M 2 (3.16) 
1-y y 2nl-y 

and 

Bj = = -£LJ*L 7 (3.17) 

3 2ti 1-y 2-Kl-y K ' 
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and Jj^ is given by 



J: 



i'j's'u 



I'J'S' 



(3.18) 



i-j/ 



I'J'S' 



The wave functions that diagonalize the left-hand side of ( 13. 15ft . and mix the physical (i = 0) 
and PV (i = 1) fermion flavors, are 



In terms of these functions, the eigenvalue problem ( 13 . 151) can be written as 



„ J i js ' 



(3.19) 



(3.20) 



Here the contribution of the zero-photon and two-photon kernels, is implicitly a func- 
tional of these new wave functions. The factors of a that appear in Ay and Bj are assigned 
the physical value and not treated as eigenvalues. The original wave functions are recovered 

aS ~ ± ■ ~ ± 

n ti± _ Ai-i,jsfij S + (—lyBjf^js 

Self-energy contributions appear in the denominators of the wave functions. 

To express the eigenvalue problem explicitly in terms of the f^ s , we first write the 
definition (13TT8D of Jj£ in a simpler form 



= I d y' d i'l E (-V^'^J^i'fAy, ox\ y', &)c$,tf, (3.22) 



where J^ s ^y s , = Jifj^j's' + Af^j's" Substitution of (I3.2ip then yields, in matrix form, 



I'J'S' V 



r (2)/«/ 



/ dy , da?X;(-i) i V 



u 0js,0j's' 


J 0js,lj's 
ljs,lj's 


( A lf 






-Ao? J 



(3.23) 



JOj's' 

iyv 

The sum over v can also be written in matrix form for the helicity components v — ±, (±) 
by the introduction of 

/ -1 0\ 
-10 
10 
\ 1 J 



A 



(3.24) 



so that 



v a,/3 



(3.25) 
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Finally, we define 

A X y By 

By — A 

as a tensor product of simpler matrices. The eigenvalue problem then becomes 



= (-i)'x s ( $ _2 f ) (3.26) 



J dy'dq'l J t^'j'Av^±iV'^±)Vj', a m''fi4s' = ( 3 ' 27 ) 

i'j's'api" 

This yields a as a function of mo and the PV masses. We then find mo such that, for 
chosen values of the PV masses, a takes the standard physical value e 2 /47r. The eigenproblem 
solution also yields the functions f£ g which determine the wave functions C^ s . From these 
wave functions we can compute physical quantities as expectation values with respect to the 
projection (12.51) of the eigenstate onto the physical subspace. 



C. Numerical solution 

The eigenvalue problem (13.271) is solved numerically. Here we discuss the method used. 
Additional details about quadratures can be found in Appendix |A] and convergence proper- 
ties are discussed in Appendix O 

The integral equations (13.271) for the wave functions of the electron are converted to 
a matrix eigenvalue problem by a discrete approximation to the integrals, as discussed in 
Appendix |A] These approximations involve variable transformations and Gauss-Legendre 
quadrature; the transformations are done to minimize the number of quadrature points 
required, in order to keep the matrix problem from becoming too large, and to reduce the 
infinite transverse momentum range to a finite interval. The resolution of the numerical 
approximation is measured by two parameters, K and N±, that control the number of 
quadrature points in the longitudinal and transverse directions. 

The integrals for the normalization and anomalous moment, (12.121) and (12.141) respectively, 
are also done numerically, but are summed over different quadrature points. These points 
take into account the different shape of the integrand that comes from the square of the wave 
functions. The values of the wave functions at these other points are found by cubic-spline 
interpolation [42| in the transverse direction. Regions of integration near the line of poles 
associated with the energy denominator require special treatment, if the poles exist, through 
quadrature formulas that take the poles into account explicitly. 

The renormalization requires finding the value of the bare mass that corresponds to the 
physical value of the coupling. This defines a nonlinear equation for the bare mass, which 
is solved with use of the Muller algorithm 42|. Finding the poles in the two-body wave 



function also requires solution of nonlinear equations, and again the Muller algorithm is 
used. 

The calculation of the anomalous moment requires computation of a transverse derivative 
of the wave functions. Because the quadrature points used for integration are not uniformly 
spaced, they are not convenient for estimating the derivative directly. Instead, the wave 
functions are first approximated by cubic splines; the derivatives are then obtained from the 
splines. 

To solve the eigenvalue problem, we treat the two-photon contributions explicitly, but 
still nonperturbatively, as corrections to the one-photon truncation with self-energy, solved 



17 



in [ll[|. We do this by considering the coupled system 

(M 2 - ml)z a /z = Yl [ dy'd^ZVf^f^/zo, (3.28) 



i'j's'a/3i" 



ftis I Z 



-J^^-V^Za/zo (3.29) 

a 

~2^ J dy ' dq/ * J iS7j's'Vj',^,i'i''l^ s >/z , 

i'j's'a/3i" 

which can be obtained from (13.11) and (I3.27p . with use of the factorization (13.121) for 
and the connection (I3.2ip between the original and flavor-mixed two-body wave functions. 

The solution is found by iteration. When the index a in f!3.28|) is equal to zero, we obtain 
an equation for m , 



m = + 



\ 



M2 "\/! E / Wtgvfl&nr^efPr,/** (3.30) 

' i'j's'a/3i" ^ 

and when a is equal to 1, we obtain an equation for z 1( 

Zl ^ Z ° = M 2 -m 2 \/f / d y' d Q± V $u>Vj',ai3,iH'JI!^ s/ /zo. (3.31) 

1 ' i'j's'afii" J 




These provide the updates of mo and Zi/zo, and (I3.29f) is solved by Jacobi iteration (42] of 
the linear system that comes from the discretization of the rearrangement 

+ ^ / dy'dq'l E J SSv%M««^Ao (3.32) 

i'j's'afii" 

a 

Only a few Jacobi iterations are performed per update of mo and ^i/^o! further inner iteration 
is unnecessary, due to the subsequent changes in m and Z\f z$. The outer iterations of the 
full system of equations is terminated when the changes in m , Zi/z , and the two-body 
wave function are all of order 10~ 6 or less. The bare amplitude Zq is obtained at the end 
by normalization. The coupling a is held fixed at the physical value; hence, this iterative 
method yields not only the two-body wave functions and one-body amplitudes but also the 
bare mass, mo- 

The eigenvalue problem must first be solved for M = 0, with the coupling strength 
parameter £ 2 adjusted to yield m = 0. This determines the value of £ 2 that restores the 
chiral limit nonperturbatively. The eigenvalue problem can then be solved for M = m e , the 
physical mass of the electron, and the anomalous moment calculated. 

The number of wave-function updates in the Jacobi iteration is small enough that the 
matrix representing the discretization of the integral equations can be computed at each 
iteration without making the calculation time too large. Thus, the matrix need not be stored, 
which allows much higher resolutions. We find that reasonable results are not obtained until 
the longitudinal resolution K is at least 50 and that there is still considerable sensitivity to 
the longitudinal resolution. There is much less sensitivity to the transverse resolution, for 
which 7Vj_ = 20 is found sufficient. 
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IV. RESULTS 



Our results for particular values of longitudinal resolution K are plotted in Fig. [3j Several 
different values are considered for the PV photon mass fit, with /x 2 fixed as a/2/ii and the 
PV electron mass ni\ set to 2 • 10 4 m e . The results are sensitive to K even for these higher 
resolutions, with greater sensitivity for the larger fi\ values. In fact, beyond /xi ~ 300m e , 
convergence is difficult to obtain for any value of K. 

The choice of value for fi2 was studied in [9]. There is no particular sensitivity to the 
choice, provided \ii is greater than \i\ and much less than m\. If Hi is less than /ii, the 
assignment of negative and positive metrics of the two PV photons must be reversed. If \ii 
is too close to mj, observables can have a strong dependence on the PV masses. 

We extrapolate the results for the anomalous moment with linear fits in 1/K. The 
estimated 5% error in the individual values, discussed in Appendix (Cj does not justify a 
higher-order fit. Given the nature of the fits, we estimate an error of 10% in the extrapolated 
values. 

The results of the extrapolations are plotted in Fig. HI Each value is close to the Schwinger 
result and independent of /ii, to within numerical error. However, there is clearly a system- 
atic tendency to be below the Schwinger result by approximately 10% as the PV photon 
mass Hi is increased. We expect that this discrepancy is caused by the absence of two poten- 
tially important contributions, the electron-positron loop and the three-photon self-energy. 
The loop contributes in perturbation theory at the same order in a as the one-electron/two- 
photon Fock states considered here and corresponds to the addition of two-electron/one- 
positron Fock states to our truncation. The self-energy contribution is higher order in a 
but earlier calculations [ll[ have shown that the two-photon self-energy is an important 
correction to the one-photon truncation. This can be seen in Fig. HJ where we reproduce 
these results for comparison. 

Figure H] also includes results obtained for the two-photon truncation when only the one- 
loop chiral constraint is satisfied. Without the full nonperturbative constraint, the results 
are very sensitive to the PV photon mass /ii. This behavior repeats the pattern observed in 
[9J] for a one-photon truncation without the corresponding one- loop constraint. The resulting 
/ii dependence is illustrated in Fig. 2 of [9J. Thus, a successful calculation requires that the 
symmetry of the chiral limit be maintained. 

V. SUMMARY 

The results of the calculation are shown in Fig. H] and compared with three other calcu- 
lations: the one-photon truncation [9j , the one-photon case with the two-photon self-energy 
contribution [llj], and the two-photon truncation with only the one-loop chiral constraint. 
The two-photon results with the correct chiral constraint are consistent with the Schwinger 
result, and therefore with experiment, to within the estimated numerical error of 10%. The 
systematic deviation below the Schwinger result is expected to be due to the absence of the 
two-electron/one-positron Fock sector and the three-photon self-energy contributions. As 
is well known from perturbation theory, cancellations exist between different types of con- 
tributions, such as between photon loops and electron-positron loops, and, therefore, it is 
not surprising for the present two-photon calculation, which does not also include electron- 
positron loops, to have a somewhat worse result than the one-photon calculation with just 
the two-photon self-energy contribution. 
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1/K 

FIG. 3. Dependence on longitudinal resolution K of the anomalous moment a e of the electron 
in units of the Schwinger term (a/27r) for the two-photon truncation. The PV masses are m\ = 
2 • 10 4 m e , [i\ = 100m e to 330m e , and \i2 = V^^i- The transverse resolution is iV|_ = 20. The lines 
are linear fits. The errors in the individual points are estimated to be 5%. 

Inclusion of an electron-positron pair in the basis is also important for the understanding 
of current covariance and of nonperturbative renormalization of the charge and photon mass. 
Future work along these lines will be to include these additional contributions. Of course, 
the calculation of the anomalous moment is not the important objective; instead, we are 
interested in testing on QED a nonperturbative method that could be applicable to QCD, to 
see how the various truncations affect the result and to be able to compare with perturbation 
theory, as a check. 

Further tests of the method within QED could include application to the calculation 
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FIG. 4. The anomalous moment of the electron in units of the Schwinger term (a/2ir) plotted 
versus the PV photon mass, fix, with the second PV photon mass, fi2, set to \f2[i\ and the PV 
electron mass mi equal to 2 • 10 4 m e . The solid squares are the result of the full two-photon 
truncation with the correct, nonperturbative chiral constraint. The open squares come from use of 
a perturbative, one-loop constraint. Results for the one-photon truncation 0] (solid line) and the 
one-photon truncation with the two-photon self-energy contribution ll|] (filled circles) are included 
for comparison. The resolutions used for the two-photon results are K = 50 to 150, combined with 
extrapolation to K = oo, and N± = 20. 



of a true bound state, positronium [43j] . Just as for the anomalous moment calculation, 
the positronium results will not be competitive with high-order perturbation theory. The 
numerical errors are large compared to the tiny perturbative corrections in such a weakly 
coupled theory. This in turn suggests another interesting test that could be done, a cal- 
culation of the anomalous moment when a is much larger that its natural value, yet small 
enough for perturbation theory to still provide a check. 

In a strongly coupled theory, such as QCD, the method may be more quantitative. For 
QCD, the PV-regulated formulation by Paston et al. |44j could be a starting point. The 
analog of the dressed-electron problem does not exist, of course, and the minimum trunca- 
tion that would include non-Abelian effects would be to include at least two gluons. The 
smallest calculation would then be in the glueball sector. In the meson sector, the minimum 
truncation would be a quark- ant iquark pair plus two gluons, which as a four-body problem 
would require discretization techniques beyond what are discussed here, since the coupled 
integral equations for the wave functions cannot be analytically reduced to a single Fock 
sector. Instead, one would discretize the coupled integral equations directly, in analogy with 
the original method of DLCQ and diagonalize a very large but very sparse matrix. 
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As an intermediate step, one can select a less ambitious yet very interesting challenge of 
modeling the meson sector with effective interactions, particularly with an interaction to 
break chiral symmetry [45j|. 
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Appendix A: Discretizations and quadratures 

The integral equations involve integration over the longitudinal momentum fraction and 
the square of the transverse momentum. The normalization, anomalous moment, and self- 
energy contributions also require integrals of this form. In each case there can be a line of 
poles gp 0le in the integrand, from the denominators of wave functions in (I3.2ip , for a range 
of values of the longitudinal momentum fraction y. The location of the line is determined 
by the energy denominator that appears in each integrand. For simple poles, the transverse 
momentum integral is defined as the principal value. For those values of longitudinal mo- 
mentum y for which the pole exists, the q\ integration is subdivided into two parts, one 
from zero to 2q^ ole and the other from there to infinity. If the pole does not exist, transverse 
integration is not subdivided. When self-energy effects are included, the location of the pole, 
if it still exists, must be found by solving a nonlinear equation numerically. We do this with 



the Miiller algorithm [42 



For the interval that contains a simple pole, the integral is approximated by an open 
Newton-Cotes formula that uses a few equally-spaced points placed symmetrically about 
the pole at qf = (2i — l)q 2 ole /N with i = 1, . . . , N and N even. This particular Newton- 
Cotes formula uses a rectangular approximation to the integrand, with the height equal to 
the integrand value at the midpoint of an interval of width 2q 2 ole /N. An integral is then 
approximated by 

f 2q l^ 2q 2 , A 

Jo iV i=i 

The equally spaced points provide an approximation to the principal value. 

This form avoids use of q\ = as a quadrature point. Such a choice is important 
for evaluating terms with two-photon kernels, where there is another pole associated with 
the three-particle energy denominator. By keeping q\ nonzero, this pole can be handled 
analytically as a principal value in the angular integration, as discussed in Appendix [Bl 

For the infinite intervals, q\ is mapped to a new variable v by the transformation [s| 

2 2 1 - (g/q ^ 

{IP /a" 



<a=^ :~z^' , (A2) 



with v in the range to 1. (If the pole exists, this transformation is shifted by 2gp 0le .) 
The PV contributions make the integrals finite; therefore, no transverse cutoff is needed. 
Only the positive Gauss-Legendre quadrature points of an even order 2N± are used for v 
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between -1 and 1, so that v = 0, and therefore q' ± = (or 2gp 0le ), is never a quadrature 
point. The points in the negative half of the range, which would be used for representing 
q\ G [— oo, 0], are discarded. One could map q\ G [0, oo] to [—1, 1] and not discard any part 
of the Gauss-Legendre range; however, the quadrature would then place points focused on 
some finite q\ value, rather than on the natural integrand peak at q\ = 0. The total number 
of quadrature points in the transverse direction is N± + N, with N = when there is no 
pole and N± typically of order 20. 

This transformation was used in [8] and was selected to obtain an exact result for the 
integral j[l/(a 2 + q 2 ) — l/(b 2 + q 2 )]dq 2 . In the present work, the scales a 2 and b 2 are 
chosen to be the smallest and largest scales in the problem, i.e. a 2 = \q 2 ole \ and b 2 = 
m\y + n\ (1 — y) — M 2 y{l — y). Here gp 0le is the location of the root of the nonlinear equation 
for the pole. If gp 0lc is negative, a pole does not exist; however, |<j^ ole | is still a natural scale 
for the integrand. 

For the normalization and anomalous moment integrals, the transverse quadrature scheme 
is based on a different transformation 



?! = a 2 ^, (A3) 
where, again, if the pole exists, the transformation is shifted by 2q 2 lc . Cubic-spline interpo- 



lation [42] is then used to compute the values of the wave functions at the new quadrature 
points. This transformation is selected to yield an exact result for the integral of l/(a 2 + q 2 ) 2 , 
which is the form of the dominant contribution to the normalization and anomalous moment. 

The longitudinal integration is subdivided into three parts when the line of poles is 
present. Two parts are symmetrically placed about the logarithmic singularity at y po i e that 
arises where the line of poles reaches q\ = 0. When self-energy effects are not included 
in the energy denominator, this occurs at y po \ c = 1 — m^/M 2 ; when self-energy effects are 
present, the location must be found by solving a nonlinear equation. The third part of the 
integration covers the remainder of the unit interval. Specifically, these intervals are [0, y po \e], 
[Vpoie, 2?/poie], and [2y pole ,l]. This structure is designed to maintain a left-right symmetry 
around the logarithmic singularity, because in the normalization and anomalous moment 
integrals (which use the same longitudinal quadrature points) the singularity becomes a 
simple pole defined by a principal-value prescription. The left-right symmetry then assures 
the necessary cancellations from opposite sides of the pole. When no pole is present, the 
longitudinal integration is not subdivided. 

The intervals are each mapped linearly to y G [0, 1] and then altered by the transforma- 
tions @ 

y{t) = t 3 {l + dt)/[l + d-(3 + Ad)t + (3 + 6d)t 2 - 4dt 3 + 2dt 4 } (A4) 

and 

t( u ) = (u + l)/2. (A5) 

The new variable u ranges between -1 and 1, and standard Gauss-Legendre quadrature is 
applied. The transformation from y to t is constructed to concentrate many points near the 
end-points of each interval, where integrands are rapidly varying. The parameter d is chosen 
such that y ~ O.Olt 3 for small t. The transformation was found empirically [8], beginning 
with a transformation constructed to compute the integral f Q [ln(y + e ) — ln(y + ei)]dy 
exactly, with eo and e\ small. The symmetry with respect to the replacements t — > (1 — t) 
and y — > (1 — y) is not necessary but is the simplest choice for restricting the coefficients in 
the denominator of (1A4D. 
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The need for a concentration of longitudinal quadrature points near and 1 is particularly 
true for the integral J, defined in (1C2j) . Although this integral can be done analytically for 
the case of the one- photon truncation discussed in [3], the integral is only implicit in the 
integral equations for the two-body wave functions discussed in Sec. [Till and must therefore 
be well represented by any discretization of the integral equations. After the transverse 
integration is performed, the integrand is sharply peaked near y = and y — 1, at distances 
of order mo/wii ~ 1CT 10 from these end-points, and needs to be sampled on both sides of 
the peaks. 

The number of points in each of the three intervals is denoted by K, which becomes the 
measure of the resolution analogous to the harmonic resolution of DLCQ [l2l ]. Thus the 
total number of quadrature points in the longitudinal direction is 3K, with K typically of 
order 50 or higher. 

For those longitudinal integrals with an upper limit less than 1, the integrand is trans- 
formed as above and given a value of zero for the points beyond the original integration 
range. 

For the normalization and anomalous moment integrals, the pole in the transverse integral 
(when it exists) is a double pole, defined by the limit j3] 

f(y,ql) 



dy dq ± 



lim -e 

e^O 2 



[m?y + -y)- M*y(l - y) + q\f 



dy I dq 2 ± f(y,k 2 ± ) 



[m 2 y + ^(l-y)-M 2 y(l-y) 



[m 2 y + -y)- M 2 y(l - y) + q\ + e] 



(A6) 



The simple poles that remain are prescribed as principal values. Of course, the limit must 
be taken after the integral is performed. 

This limiting process is taken into account numerically by using a quadrature formula 
that is specific to this double-pole form. On the interval [0,2gp Ole ], the quadrature points 
are chosen to be the same as those used for the integral equations, which are qf = (2i — 
1, . . . , N, as given above. The interval is divided into N/2 sub intervals 

0, 1, . . . , (N — 2)/2, each containing two of the quadrature 



l)<Zpole/^ With i 



r 4m „2 4(m+l) 2 1 
I JV Vpole' N Vile.!' wlul m 



points. The quadrature formula for such a subinterval is taken to be 



4(m + l) 2 
N "pole 



4m 2 
N ^polc 



dq\ 



W 2m +lf( 



(Am + 1) 
N 



2 \ r/(4wl+3) o / A jt\ 

tfpole) + »2m+2/( ^ g po l e ), (A7) 



where the integral on the left is defined by the limit formula in ( ]A6j) when the pole is in the 
subinterval. The weights Wi are chosen to make the formula exact for / = 1 and / = q\ 
on each individual q\ subinterval. For these numerator functions, the limit in ( 1A6[) can be 
taken explicitly. The weights are then found to be u>at/2 = wnI2+\ = — ^/2?p i e , for the 
quadrature points on either side of the pole. For all other points, the weights are given by 
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The integral from to 2q^ ole is obtained by summing over the individual subintervals. 

For the self-energy contribution (13.101) . which is expressed in terms of the integrals Iq, 
Ji, and J = M 2 Iq given in (ICip and (IC2I) . the transverse integral is done analytically. Only 
the longitudinal integral is done numerically, by the scheme discussed above with resolution 
K = 30. 



Appendix B: Angular integrals 

Calculation of the two-photon kernels requires the integrals 

r r d(p e " m0 ran 

n ~J 2vrD + Fcos0' 1 } 

first defined in Eq. (13. 7p . with D and F given in ( I3.8p . Here the original integration variable 
0' has been shifted by the independent angle 0, and the prime then dropped for simplicity 
of notation in this Appendix. The factor F is always positive, but D can be negative. If 
the bare fermion mass mo is less than the physical mass m e , we can have \D\ < F; in this 
case, X n is defined by a principal value, as in the one-photon sector. If either photon has 
zero transverse momentum, F will be zero, and any pole due to a zero in D will not involve 
the angular integration. The numerical quadrature is chosen to never use grid points where 
a photon transverse momentum is zero, so that the principal-value prescription can always 
be invoked for the angular integral, where it is easily handled analytically. 

The imaginary part of X n is zero. This follows from the even parity of the denominator 
and the odd parity of sinn0. As a consequence, X_ n = Z n , and we evaluate (IB1[) for only 
nonnegative n. 

The real part is nonzero and most easily calculated from combinations of the related 
integrals 

X r d4> C ° Sn0 (B2) 
n J 2tt D + F cos cf>' 1 ' 

Of course, for n = and 1, the two integrals are identical. For n = 2 and 3 we have 
cos 20 = 2 cos 2 — 1 and cos 30 = 4 cos 3 — 3 cos 0. Therefore, the integral combinations 
are 

= 2^o, Xi = Xi, X 2 = 2X 2 — X , X 3 = 4X 3 — 3Xi. (B3) 

Larger values of n do not appear in the two-photon kernels. 

The integrals X n are connected by a simple recursion for n > 0: 

= f* d0cos" -1 </>(D + F cos </> - D) = 1 /* # D 

J 2tt F D + F cos F J 2tt v F y 1 

The first term is zero when n is even. For n = 1, it is 1/F, and for n = 3, this term is 1/2F. 
The only other integral that must be evaluated directly is X =X . 

The determination of Xq, with or without the presence of poles, is conveniently done by 
contour integration around the unit circle in terms of a complex variable z = e 1 ^. We then 
have 

If dz 1 / dz 

G = mFj zi + 2§z + l = ™Fj (z-z + )(z-z.Y (B5) 
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There are simple poles at 



When D is greater than F, one pole, z + , is inside the contour and the other outside, as 
illustrated in Fig. [5j Evaluation of 2iri times the residue yields 

y 




FIG. 5. Integration contour for evaluation of Xq. The locations of the poles at z± depend upon the 
magnitude and sign of D/F. The semicircles are used when \D\/F < 1 and the poles are on the 
contour. 



X = ; for D > F. (B7) 

VD 2 - F 2 V ; 

Similarly, when D is less than —F, the pole at z_ is inside the contour, and we have 

X = ; 1 for D < -F. (B8) 

VD 2 - F 2 V ; 

When \D\ is less than F, the poles move to the contour, and the integral is defined by the 
principal value. This is evaluated by distorting the contour to include semicircles of radius e 
around each pole, as shown in Fig. G3, and subtracting the contributions from the semicircles 
after taking the e — > limit. The choice of inward semicircles makes the integral around 
the closed contour simply zero. For the semicircle around z±, we have z — z± + ee ld and a 
contribution, as e goes to zero, of 

2 f eie iS dd ^ ±JdB 



iF J ee ie {ee ie =f 1iyJ\ - D 2 / F 2 ) Fy/l - D 2 /F 2 ' 
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Thus the contributions from the two semicircles are of opposite sign and cancel, so that the 
net result is also zero. Therefore, we have 



i 



X 



VD 2 -F 2: 

o, 

1 



D > F 

\D\ < F 
D < -F. 



(BIO) 



The case where D equals F represents an integrable singularity for the transverse momentum 
integrations and can be ignored. When F is zero, we have simply 
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(Bll) 



When F/D is small, the expressions for the integrals X n are best evaluated from ex- 
pansions in powers of F/D, to avoid round-off errors due to cancellations between large 
contributions. The expansions used are 
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Appendix C: Numerical convergence 

The primary constraint on numerical accuracy is the error in the estimation of the inte- 
grals in the integral equations for the wave functions and in the expressions for the normal- 
ization and the anomalous moment. This accuracy is determined by the choice of quadrature 
scheme, discussed in Appendix [A] and the resolution, controlled by the longitudinal param- 
eter K and transverse parameter N±. The other numerical parts of the calculation are 
iterated to what is effectively exact convergence, with remaining uncertainties much smaller 
than the errors in the numerical quadratures. 

For the one-photon truncation, discussed in [3] and all the integrals can be done 
analytically. This makes the one-photon problem a convenient first test for numerical con- 
vergence. The key integrals are Iq, Ji, and J, defined as 

U M>) = / *tM V (Zip . , (C1) 



j(m 2 ) - I v (-i)^ ™* + kj 

3 l-v y 
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TABLE I. Dependence on longitudinal resolution K of the integrals Iq, I\, and J, denned in (|Clj) 
and (|C2p of the text and computed according to the quadrature scheme described in Appendix El 
The bare-electron mass is mo = 0.98m e . The PV masses are m\ = 2 • 10 4 m e , fi\ = 200m e , and 
fj,2 = V^fJ-i- The transverse resolution is iVj_ = 40. 



K 


h{m 2 e ) 


Ii{m 2 e )/m e 


J (ml) 1 ml 


5 


-7.113 


-13.530 


-11606. 


10 


-6.2586 


-10.6182 


932.3 


15 


-6.2641 


-10.7354 


-26.487 


20 


-6.2645 


-10.7327 


-6.7126 


25 


-6.2645 


-10.7328 


-6.3982 


30 


-6.2645 


-10.7328 


-6.4401 


exact 


-6.2645 


-10.7328 


-6.2645 








10 



20 



30 



K 



FIG. 6. Dependence on longitudinal resolution K of the integrals lo, I\, and J, defined in (|C1|) 
and (|C2p of the text and computed according to the quadrature scheme described in Appendix lAl 
The values plotted are ratios to the exact values, listed in Table HI The bare-electron mass is 
mo = 0.98m e . The PV masses are m\ = 2 • 10 4 m e , fii = 200m e , and ^2 = V^fJ-i- The transverse 
resolution is N± = 40. 



Tables [I] and [TT] and Figs. [6] and [7] summarize results for numerical calculation of these 
integrals. They show that Jo and I± are well approximated for a wide range of resolutions, 
but J is particularly sensitive to the longitudinal resolution K and requires that both K 
and iVj_ be on the order of 20 or larger. At these resolutions, J is approximated with an 
accuracy of about 4%, and this then becomes a minimal estimate of the accuracy of any of 
the results. 
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TABLE II. Same as Tabled! but for the dependence on transverse resolution N±. The longitudinal 
resolution is K = 30. 



N± 


h{ml) h(m 2 e )/m e 


J{m 2 e )/m 2 e 


10 


-6.2646 


-10.7324 


-2.9845 


15 


-6.2645 


-10.7327 


-7.5428 


20 


-6.2645 


-10.7327 


-6.9137 


25 


-6.2645 


-10.7328 


-6.6961 


30 


-6.2645 


-10.7328 


-6.5712 


35 


-6.2645 


-10.7328 


-6.4922 


40 


-6.2645 


-10.7328 


-6.4401 


exact 


-6.2645 


-10.7328 


-6.2645 
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FIG. 7. Same as Fig. [HJ but for the dependence on transverse resolution N±. The longitudinal 
resolution is K = 30. For these resolutions, the values for Jo and I\ are nearly exact, and the 
plotted points for the ratios to the exact values are at the same places; only J shows variation. 



As expected, the results for the one-photon truncation, if computed numerically, converge 
to better than 1% at the same resolution, of K = 20 and N± = 20, as can be seen in Tables ILTT1 
and IIVI and Figs. [S] and However, the results with the self-energy contribution, shown 
in the same tables and figures, require K ~ 25 before nearing convergence. Although the 
exact answer is not known in this case, K = 20 is clearly insufficient, but K ^> 25 yields a 
reasonable result with an error on the order of 1%. 

The two-photon truncation incorporates numerical approximations to the integrals Iq, I\, 
and J through the action of the zero-photon kernel J^°\ in Eq. (13.91) . and approximations 
to the self-energy contribution, also in Eq. (13.91) . Thus, the minimum resolution for the 
two-photon calculation would appear to be approximately K = 25 and N± = 20; however, 
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TABLE III. Dependence on longitudinal resolution K of the bare mass mo and anomalous moment 
a e of the electron in units of the physical mass m e and the Schwinger term (a/27r), respectively, 
for both the one-photon truncation, when solved numerically, and the case with the self-energy 
included. The value of the bare mass is obtained as the value that yields a physical value of the 
coupling constant a; for K = 15, two solutions are found. The PV masses are m\ = 2 • 10 4 m e , 
Hi = 200m e , and fX2 = V^i- The transverse resolution is iV_|_ = 40. 





one-photon 


with self-energy 


K 


mQ/m e 2ira e /a 


mo/m e 2ira e /a 


5 


4.4849 0.09921 


4.4845 0.09444 


10 


1.7445 0.22992 


1.7106 0.22620 


15 


1.07487 0.64482 


1.07481 0.59092 


15 


0.98906 1.12763 


0.97537 1.07444 


20 


0.98223 1.15430 


0.99028 0.90337 


25 


0.98240 1.15382 


0.98295 0.99242 


30 


0.98241 1.15567 


0.98245 1.00612 



1.2 
1.0 

0.8 - 

0.4 
0.2 
0.0 







T 



• one-photon truncation 
▼ with self-energy 



10 
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FIG. 8. Dependence on longitudinal resolution K of the anomalous moment a e of the electron in 
units of the Schwinger term (a/2ir) for both the one-photon truncation, when solved numerically, 
and the case with the self-energy included. The PV masses are mi = 2 • 10 4 m e , fii = 200m e , and 
li2 = V2/J-1- The transverse resolution is iV_|_ = 40. 



we find that K must be at least 50. We extrapolate from the one-photon and self-energy 
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TABLE IV. Same as Table HIH but for the dependence on transverse resolution N±. The longitu- 
dinal resolution is K = 30. 



N± 


one-photon 


with self-energy 


mo/m e 2i:a e /a 


mo/m e 2ira e /a 


10 


0.98023 1.16121 


0.97918 


15 


0.98297 1.15317 


0.98307 0.99071 


20 


0.98268 1.15380 


0.98275 0.99776 


25 


0.98256 1.15484 


0.98261 1.00130 


30 


0.98248 1.15399 


0.98253 1.00348 


35 


0.98244 1.15404 


0.98248 1.00501 


40 


0.98241 1.15567 


0.98245 1.00612 
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1.1 



• ••••• 



5 i.o 

CM 
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FIG. 9. Same as Fig. [HI but for the dependence on transverse resolution N±. The longitudinal 
resolution is K = 30. 



calculations to estimate an error of 5-10% for the two-photon truncation. 
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